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Transverse and parallel static susceptibilities of in-plane uniaxial anisotropic ferromagnetic films 
are calculated using a Heisenberg model within the framework of many-body Green's function 
theory. The importance of collective magnetic excitations, in particular in the paramagnetic regime, 
is demonstrated by comparing with mean field calculations. The paper extends the work of Jensen 
et ali on the monolayer with spin 1/2 to the multilayer case with arbitrary spin. 

1. Introduction 

In a recent paper, Jensen et ali reported on the measurement of the magnetic susceptibility of a bilayer Co film 
with an in-plane uniaxial anisotropy. On the basis of a Heisenberg model they determined the isotropic exchange 
interaction and the magnetic anisotropy within the framework of a Green's function theory by fitting the interaction 
parameters of the theory to the measured susceptibilities along the easy and hard axes in the paramagnetic regime, 
assuming a spin S — 1/2. In this paper, we generalize the theoretical treatment to the multilayer case and to spins 
S>l/2. We organize the paper as follows. In Section 2 we explain the model and the Green's function formalism for 
its solution. Section 3 displays the numerical results. In the final Section 4 we summarize the results and present our 
conclusions. 



2. The model and the Green's function formalism 

We consider a Hamiltonian consisting of a ferromagnetic isotropic Heisenberg exchange interaction with strength 
( Jki > 0) between nearest neighbour lattice sites, a uniaxial in-plane exchange anisotropy in the z-direction with 
strength (D k i > 0), and an external magnetic field B = (B x ,0, B z ) confined to the film plane: 

W = ~ o ^ J ki(S^S^ + SISi) - - ^ D kl SlSi 

<kl> <kl> 

Y,( BX \( S k+ S k)+B z S* k ). (1) 

Here the notation S k = S k ± iS\ is introduced, where k and I are lattice site indices and < kl > indicates summation 
over nearest neighbours only. In keeping with reference^ we do not consider the dipole coupling, since it is almost 
isotropic for the in-plane situation. We note, however, that the formalism is capable of handling this coupling, 
if so desired, and we refer the reader to references^, where we treated the reorientation of the magnetization of 
ferromagnetic films with anisotropies normal to the film plane. Also in keeping with reference^, we choose the 
exchange anisotropy over the single-ion anisotropy, which has the advantage of being simpler to handle in the Green's 
function theory. Moreover, we have shown in4 that once the strength of the exchange anisotropy is fitted appropriately, 
the magnetization as a function of the temperature and film thickness behaves very similar to that calculated from 
the single-ion anisotropy, which may appear somewhat surprising, since the anisotropies originate from very different 
physical mechanisms. In this paper we shall restrict ourselves to a simple cubic lattice. 

In order to generalize the treatment of - to general spin S, we need the following Green's functions 

Gt^ n (u J ) = ((s^(s z r(s-r))^ , (2) 

where a = (+, — , z) takes care of all directions in space, r\ = ±1 refers to the anticommutator or commutator Green's 
functions, respectively, and n > l,m > are positive integers, necessary for dealing with higher spin values S. For 
n = 1 and m = we will recover the equations oft. 

The exact equations of motion for the Green's functions are 

wG&™» = + (([S?,HU (S*r i (Srr))^ n (3) 
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with the inhomogeneities 

^;7H[sf,(^ns7m (4) 

where (...) = Tr(...e~P n )/Tr(e~P' H ) denotes the thermodynamic expectation value. 

After solving these equations, the components of the magnetization can be determined from the Green's functions 
via the spectral theorem. A closed system of equations is achieved by decoupling the higher-order Green's functions 
on the right hand sides. For the exchange-interaction and exchange-anisotropy terms, we use a generalized Tyablikov- 
( RPA-) decoupling 

(S*r(S-) n )) v c (ST)G^ + (S%)G^ n . (5) 

We do not try to do better than RPA because we have shown in earlier work^, by comparing RPA with "exact" 
Quantum Monte Carlo calculations that the former is quite a good approximation in simple cases. 

We now proceed with the formulation of the theory for the multilayer case. After a Fourier transform to momentum 
space, we obtain, for a film with N layers, 3N equations of motion for a 37V-dimensional Green's function vector G" m : 

(ujI - T)G™" = A™", (6) 

where 1 is the 3iV x 3-/V unit matrix. The Green's function vectors and inhomogeneity vectors each consist of N 
three-dimensional subvectors which are characterized by the layer indices i and j 



G™(k, W ) 
G£»(k, W ) = ( G£T(k, W ) 



a mn 
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mn 
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(7) 
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The equations of motion are then expressed in terms of these layer vectors and the 3x3 submatrices of the 
3N x 3N matrix T 
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(8) 



When performing the decouplings according to equation (0 , the T- matrix reduces to a band matrix with zeros in the 
Tij sub-matrices, when j > i + 1 and j < i — 1. The diagonal sub-matrices Tu are of size 3x3 and have the form 





1 TJX 





- H J 

1 TT1 

2 H i 







(9) 



where 



HI = Z l + (S*)Mq~Jk) , 
Z t = B z +D ii q(S?) + (J i , i+ i+D i , i+1 )(Sf +1 ) 

+(J M _ 1 +A,i-i)(S?_ 1 ) , 
Hf = B a: + (5f>J»( g -7k)+J M+ i(S? +1 > + J i , i _ 1 (5f_ 1 ) 
Hf = Hf - (Sf)D iak . 



(10) 



For a square lattice and a lattice constant taken to be unity, 7k = 2(cos/cj; + cosfc z ), and q = 4 is the number of 
intra-layer nearest neighbours. The mean field (MFT) results, which we use later for comparison with Green's function 
theory, are obtained by putting 7^ = 0; i.e. only the number of nearest neighbour counts, whereas RPA introduces a 
momentum dependence on the lattice under consideration. 
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Note that because the momentum dependence in Hf stems from the exchange anisotropy, Hf ^ Hf , which prevents 
a naive extension of the formalism of reference^. 

The 3x3 off-diagonal sub-matrices T^- for j = i ± 1 are of the form 

/ -Jij(Sf) (Jij + Dij)(Sf) \ 

r«= o j^Sf) -(j ij + Dij) {Sf) . (ii) 
V \Jn(Sf) -\Ja(s?) o J 

The treatment of multilayers is only practicable when using the eigenvector method developed in reference^. The 
essential features are as follows. One starts with a transformation, which diagonalizes the T-matrix of equation JBJ 

LrR=f2, (12) 

where 1~2 is a diagonal matrix with eigenvalues lo t (t = 1, 3iV), and the transformation matrix R and its inverse 
R 1 = L are obtained from the right eigenvectors of T as columns and from the left eigenvectors as rows, respectively. 
These matrices are normalized to unity: RL=LR=1. 

Multiplying the equation of motion JBJ from the left by L and inserting 1=RL one finds 

(cjI - ft)LG m " = LA™. (13) 

Defining Q™ n = LG mn and A™ n = LA™" one obtains 

{uji-n)g™ n = a™ 1 . (14) 

G™ n is a vector of new Green's functions, each component r of which has but a single pole 

jmn,r 

gmn,r = _n _ (15) 

' ui — Ul r 

This is the important point because it allows application of the spectral theorem, e.gi, to each component separately. 
We obtain for the component r of correlation vector C mn = LC mn ( where C mn = {(S z ) m (S-) n S a ) with (a = 

+.-*)) 

Qmn.j _ _2H + i ( l_ r?) il imtJ ^!E±l^ (16) 

We emphasize that when (j] = —1), the second term of this equation, which is due to the anticommutator Green's 
function, has to be taken into account. This term occurs for uj t = and can be simplified by using the relation 
between anticommutator and commutator 

Aj=+1 = - 4 7 7 =-l + 2C = L (A ?J= _ 1 + ), (17) 

where the index r = refers to the eigenvector with lu t = 0. 

The term L,qA™ 1 _ 1 — vanishes due to the fact that the commutator Green's function has to be regular at the 
origin 

lim LUG*'™ = 0, (18) 

which leads to the regularity conditions: 

H x A+^l + H x A-ft\ + 2H z A z ^ m 1 = 0. (19) 

For details, see reference^. 
This is equivalent to 

Lo^-i - 0, (20) 
because the left eigenvector of the T-matrix with eigenvector zero has the structure (see also equation l|27|l below) 

L oc (H X ,H X ,2H Z ). (21) 
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For more details concerning the use of the regularity conditions, see refsm^. 

The equations for the correlations are then obtained by multiplying equation l|l(jfl from the left with R and using 
equation JSJJJl; i.e. 

C = R£LA + R L C, (22) 

where £ is a diagonal matrix with matrix elements £{j = 5ij(e^ 3uJi — 1) _1 for eigenvalues 7^ 0, and for eigenvalues 
Ui = 0. 

A problem associated with this equation is that the exchange anisotropy introduces a momentum dependence into 
the projection operator R0L0. Consequently, when the Fourier transform to real space is performed, the projector 
cannot be taken out of the integral as is possible in the case of the Anderson Callen decoupling of the single-ion 
anisotropj^, where the projector turns out to be momentum independent. The solution is to eliminate one component 
of the projector by a transformation of equation (|22|l . which is sufficient to establish the integral equations of the 
eigenvector method. 

The adequate transformation is found to be 
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with T- X T = 1. 

Applying this transformation to equation ill' 21) considered as a monolayer problem 

T _1 C = T^RfLTT^A + T^RoLqTT^C (24) 

transforms the second component of the vector T _1 RoLoTT _1 C to zero. This can be seen when inserting the explicit 
expressions for the eigenvalues and eigenvectors for the monolayer. The eigenvalues of the T-matrix in this case are 

ui = 0; w 2 ,3 = ±£k = ±Vh z H z + H x H x , (25) 
the right eigenvectors are the columns of the matrix 



R 



/ HZ -(tk+fl') 

El (^z H *1 ~{^+h z ) ] ) (26) 
\T T 




and the left eigenvectors are the rows of the matrix 

1 



2H X H Z 2H X H Z 4H Z H Z 

(e k + H z )H x (e k -H z )H x 2H X H X I . (27) 



4e? , 

k \ (e k -H z )H x -(e k + H z )H x 2H X H X 

The second row of the transformed equation (|24|l . together with the regularity conditions (|19|) . leads to one integral 
equation for the correlations for each (m, n)-pair. 

The eigenvector method can immediately be generalized to the case of N layers by transforming equation 1211) 
(extended to 37V-dimcnsions) with a 37V x 3A-matrix T _1 having the 3x3 T~ x sub-matrices (12311 on the diagonal. 

Before showing numerical results for S > 1/2 and for multilayers we derive the equations for the monolayer and 
S = 1/2 of reference^ from equation (|24(1 . 

Using equations (|26I27JI , one obtains from the second row of equation (|24|l for the monolayer and for general spin 

\{A+^[ + A-^[)^Y, % coth (/3e k /2) (28) 



and from the first row 



9 T-f x T-f x 



H z r, — 1 H z 

A-fZ^Yl % ™th(/W2). 



2 V N ^ H z 

k 
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The equations for 5 = 1/2 are obtained from these equations for n = l,m — 0, i.e. C ,01 = (S S ) = 0, 
C+' 01 = (S-S+) = 1/2 - (S z ), A~£\ = 0, A+f\ = 2(S Z ), C z01 = (S-S z ) = 1/2(S X ), and from the regularity 
condition G3j) (S z ) = {H z / H*)(S X ): 

^H^>^£-^coth(/?e k /2), (30) 

k 

and 

^ = ^>iE|r coth (/3ek/2) ; (31) 

k 

which are the equations used in refi 

For S > 1/2 and for the multilayer case, the eigenvectors have to be calculated numerically and the system of integral 
equations obtained from the Fourier transform of equations (|24|l to real space has to be solved self-consistently, which 
is done by the curve-following method described in detail in£. 



3. Numerical results 



3.1. The monolayer with arbitrary spin S 

We start by showing results for the magnetizations (S z ) and (S x ) with respect to the easy and hard axes of 
a monolayer with various spin values S. It turns out that one obtains fairly universal curves when scaling the 
parameters of the model as J/S(S + 1) = J, D/S(S + 1) = D, and B x ( z ' / S = B x<yZ) . In the present paper we restrict 
ourselves to the case where the exchange interaction and exchange anisotropy parameters are the same for all layers 
and interlayer couplings. The program is, however, written in such a way that different parameters can be easily 
introduced. 

Using the scaled variables, the Curie temperatures Tc(S) collapse to a single universal value (the same for each 
spin S) both in mean field theory (MFT) and in the random phase approximation (RPA). 



RPA+MFT: <S Z >/S, 1 00*<S*>/(S+1 ) 
monolayer: S=1 12 ,1 ,3/2,2,3,4,6,1 3/2 




FIG. 1: The 'universal' magnetizations (S z )/S and 100 * (S X )/(S + 1) of an anisotropic ferromagnetic Heisenberg monolayer 
for a square lattice are shown as functions of the temperature for S = 1/2, 1, 3/2, 2, 3, 4, 6, 13/2. Comparison is made between 
Green's function (RPA) and mean field (MFT) calculations using the exchange interaction J = 75, the exchange anisotropy 
D = 3.75 (corresponding to J = 100, D = 5 of the S=l/2 case of Refii) and small magnetic fields B x = B z = 0.01. 



This is shown in Fig.l, where MFT and RPA results for (S z )/S and 100 x (S X )/(S + 1) are displayed as functions 
of the temperature for the spin values S = 1/2,1,3/2,2,3,4,6,13/2. For the exchange interaction and exchange 
anisotropy, we use the same parameters as in Ref. 1 , which is the S=l/2 monolayer case in our investigations. The 
calculations for Fig.l are done with small fields (B x = B z = 0.01); this stabilizes the numerical algorithm. Whereas 
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{S z )/S in RPA is universal over the whole temperature range, the corresponding MFT curves as function of the 
temperature split somewhat, reaching a saturation for large spin values S (approaching the classical limit). The 

curves 100 x (S x ) /(S + 1) have the same universal value in MFT and RPA for T < ^ PA ( MFT ) (Because of the very 
small B x held we introduced the factor 100 to make the curves visible). The reason why the values for (S X )/(S + 1) 
coincide below the Curie temperature is that the magnetization in x-direction depends only on the number of nearest 
neighbours. This can be understood from equation 1191) from which one obtains with (n = l,m = 0) 



X 1 B^oH* Dq DqS 



(32) 



This explains the universality of (S X )/(S + 1) = B x /(Dq), where q = 4 is the number of nearest neighbours for the 
square monolayer. The fact that the universal Curie temperature T PPA (S) is only about one half ofT^, IFT (S) for the 
monolayer is due to the action of collective excitations (magnons=spin waves), which are completely absent in MFT. 

Spin waves also have significant effects on the susceptibilities (in particular in the paramagnetic regime T > Tq) 
with respect to the easy (xzz) and hard (Xxx) axes. The susceptibilities are calculated as differential quotients 



X« = ((S'(B'))-(S*(0)))/B* 
(S X (B X )) - (S x {0)))/B x , 



X, 



(33) 



where the use of B z = B x = 0.01/5 turns out to be small enough to get good numerical estimates of the quotients; 
smaller fields would only be necessary to get better estimates close to the divergence of Xzz at Tc; however, the errors 
in the inverse susceptibilities x'zl and Xxx at this point arc not noticeable in the figures. The inverse susceptibilities as 
functions of the temperature can be brought into near coinicidence with a single universal curve if they are multiplied 
with a factor S(S + 1), especially in the paramagnetic regime. In Fig. 2 we compare RPA and MFT calculations for 

RPA+MFT inverse susceptibility X Z 1 Z *S(S+1 ) 
monolayer S=1 /2,1 ,3/2,2,3,4,6,1 3/2 
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FIG. 2: 'Universal' inverse susceptibilities \zz * S{S + 1) along the easy axis of an anisotropic ferromagnetic Heisenberg 
monolayer for a square lattice as functions of the temperature for S = 1/2, 1, 3/2, 2, 3.4, 6, 13/2. Comparison is made between 
Green's function (RPA) and mean field (MFT) calculations. 



the inverse susceptibility x.zz~S(S + 1) along the easy axis. As in Fig. 1, there is a shift in the Curie temperatures 
in going from RPA to MFT. For T < Tc, the MFT result behaves more universally than that of RPA; above the 
Curie temperature, both results are nearly universal. Xzz 1 ^ 1 ^^ * S(S + 1) is a straight line oc (T — T^, IFT ) like a 
Curie- Weiss law. For S = 1/2, one finds analytically from equation (|30() in the limit (S z ) — » 0,T — > large, that 

S(S + l) X ^ MF ^ = h(T-T^ T ), (34) 

where T^ IFT = J + D. The inverse RPA susceptibility x7z (RPA) * S(S + 1) is curved for T > T PPA due to magnon 
effects. This is a behaviour known from isotropic bulk ferromagnets, but the effect is significantly stronger for a 
monolayer. 
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RPA+MFT: inverse susceptibility X X 1 X *S(S+1 ) 
monolayer: S=1 /2,1 ,3/2,2,3,4,6,1 3/2 



300 



+ 200 
CO 



CO 



X 



100 



RPA 



•' MFT 



20 40 60 80 100 120 140 160 180 200 

temperature 



FIG. 3: The 'universal' inverse susceptibilities Xxx * S(S + 1) along the hard axis of an anisotropic ferromagnetic Heisenberg 
monolayer for a square lattice are shown as functions of the temperature for S = 1/2, 1, 3/2, 2, 3, 4, 6, 13/2. Comparison is made 
between the results of Green's function (RPA) and mean field (MFT) calculations. 



An analogous universality is obtained for the inverse suceptibility x X x * S(S + The results are shown in Fig. 3. 
Contrary to the curve for x7z> the hard axis susceptibility does not go to zero at T — Tc- For T < Tc one has the 
same universal constant in RPA and MFT, which can be calculated analytically for the monolayer Xxx *S(S+l) = Dq 
from equation (|32|l . The slopes of the curves for T > Tc are, however, different. The MFT again yields a straight line, 
whereas the RPA result is curved and approaches a straight line only for very large T. Hence, we see that owing to the 
scaling properties, it is not necessary to do calculations for each spin value separately. It suffices to do calculations 
for one spin value and then to apply scaling to obtain the results for other spin values. 



3.2. Multilayers at fixed spin S 

Next we discuss the multilayer case for fixed spin. We use the example of spin 5 = 1/2 (We have also considered 
multilayers with spins S > 1/2. The results scale with respect to the spin as in the monolayer case). 

Curie temperatures as function of the layer thickness are shown in Fig. 4. The difference between RPA and MFT 
is largest for the monolayer, where T§ PA ~ 0.60 T^ IFT and shrinks to T c lPA ~ 0.80 T C IFT for a film with N=19 
layers, where one is approaching the bulk value. 

To further emphasize the difference between RPA and MFT, we show in Figs. 5 and 6 the inverse susceptibilities 
Xjz an d Xxx as functions of the temperature. To avoid cluttering the figures, we plot only the results of the monolayer 
(N=l) and the film with the maximum number of layers (N=19), which is close to the bulk limit because the Curie 
temperatures saturate with increasing number of layers N. In each case, we observe the shift in the Curie temperatures 
between RPA and MFT corresponding to Fig. 4. 

Above the Curie temperature, the inverse susceptibilities are straight lines in the mean field case and curved lines 
in RPA. The slopes of the curves, however, are different for each film thickness. This is seen most clearly if one 
normalizes the temperature scale to the Curie temperatures Tq(N). The slopes in MFT increase with increasing film 
thickness, and in RPA the curvature decreases with increasing number of layers, as shown in fig. 7. 

For T < Tc the inverse susceptibility Xxx i s constant, having the same value in MFT and RPA but depending 
on the number of layers. The reason for the layer-dependence is that, with increasing film thickness, the number of 
nearest neighbours increases and therefore one has an increase in the inverse susceptibility. For the square lattice 
monolayer and bilayer, the values can be calculated analytically from the regularity condition 11911 : Xxx(N = 1) = Dq 
and Xxx(N = 2) = D(q + 1), with q — 4 for a square lattice. From the value of Xxx at T = Tc, one can obtain 
an estimate of the exchange anisotropy strength parameter D, which, together with the Curie temperature, which 
depends on the exchange interaction strength J and on the exchange anisotropy strength D, affords an estimation of 
J. As the number of layers N increases, the relative weights of the layers having two neighbouring layers increases, 
since it is only the sites in the surface layers which are restricted to having one nearest neighbour in the next layer; 
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Curie temperatures: Tg FT T£ PA 
as function of the layer thickness 

160 | . 1 . 1 . 1 . 1 . 1 . 1 . 1 . r- 




number of layers 



FIG. 4: Curie temperatures of ferromagnetic films with spin 5 = 1/2 are shown as functions of the film thickness for MFT 
(Tg IFT ) and RPA (T§ PA ). The difference T^ FT - T§ PA shrinks by a about a factor of two when going from the monolayer 
to the bulk limit. 




FIG. 5: The inverse susceptibilities \zz * S(S + 1) along the easy axisof a ferromagnetic film with spin S = 1/2 for RPA and 
MFT are shown as a function of the temperature for a monolayer (N=l) and a multilayer (N=19). 



hence \xx slowly increases (see Fig. 8). 



4. Summary and conclusion 



We have generalized the many-body Green's function treatment for calculating in-plane anisotropic static suscepti- 
bilities of ferromagnetic films to arbitrary spin S and to multilayers. In particular, we have emphasized the difference 
in the results from a Green's function theory (RPA) and mean field theory (MFT), pointing out the importance of 
spin waves, which are absent in MFT. All results discussed below refer to a simple cubic lattice. 

By introducing scaled variables, we were able to show that the magnetic properties of thin ferromagnetic films 
with in-plane anisotropy manifest nearly universal behaviour. Plotting (S z )/S and (S X )/(S + 1) as functions of 
the temperature reveals a nearly universal behaviour over the whole temperature range for RPA, whereas (S z )/S 
for MFT shows a small dependence on S. The main difference between RPA and MFT is in the universal Curie 
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inverse susceptibilities X„ 1 X *S(S+1 ) 



RPA and MFT for N=1 and N=1 9 layers spin S=1 12 




100 200 300 



temperature 

FIG. 6: The inverse susceptibilities \xx * S(S + 1) along the hard axis of a ferromagnetic film with spin S = 1/2 for RPA and 
MFT are shown as a function of the temperature for a monolayer (N=l) and a multilayer (N=19). 



inverse susceptibility X X 1 X *S(S+1 ) 



MFT+RPA for N=1 ,2,3,5,1 1 ,1 9 layers spin S=1 12 




12 3 



T/T Curie (N) 

FIG. 7: The inverse susceptibilities \xx * S(S + 1) along the hard axis of a ferromagnetic film with spin S =1/2 for RPA and 
MFT are shown as a function of the reduced temperature T /Tcuru(N) for N=l,2,3,5,ll,19 layers. 



temperature which, for the monolayer, is nearly a factor of two larger in MFT than in RPA, T^ IFT ~ 2 * Tq PA , due 
to spin wave effects. The hard-axis magnetization (S x ) has the same constant value in RPA and MFT for T < Tq 
because it depends only on the number of nearest neighbours and not on the momentum of the lattice. The inverse 
susceptibilities along the easy and hard axes also behave universally when scaled as x7z * S(S + 1) and Xxx *S(S+1), 
particularly in the paramagnetic regime. The difference between MFT and RPA consists in the shift of the Curie 
temperatures and the behaviour in the paramagnetic region (T > Tc). Whereas the MFT inverse susceptibilities are 
linear in (T — Tq ) (a Curie- Weiss-like behaviour), the RPA susceptibilities are curved, owing to spin- wave effects, 
and approach a straight line only asymptotically for very large temperatures. As long as one scales with respect to 
the spin S there is no qualitative change in the physical picture from that discussed in Refi for spin 5 = 1/2. It is 
not necessary to perform calculations for each spin value separately. Instead, it is sufficient to calculate the results 
for one spin value and to apply scaling. This is one of the main results of the paper. 

For the multilayers at a fixed spin S, the Curie temperature increases with increasing film thickness, approaching 
the bulk limit around a film consisting of N ~ 19 layers. The difference between the Curie temperatures for MFT and 
RPA decreases with increasing film thickness from T^ 1FT /T^ PA {N = 1) ~ 2 to T£ IFT ' /Tg PA {N = 1) ~ 1.3, which 
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FIG. 8: The inverse susceptibilities \xx * S(S + 1) along the hard axis at the Curie temperatures of a ferromagnetic film with 
spin S — 1/2 for RPA are shown as a function of the film thickness. 



shows that the spin wave effects are strongest for the monolayer. In MFT, the inverse susceptibilities show a linear 
Curie- Weiss behaviour for T > Tc, whereas the RPA results are curved. When plotting the inverse susceptibilities as 
a function of the normalized temperatures T/Tq(N) the slopes of the straight lines of MFT increase with increasing 
layer number, whereas the curvatures of RPA decrease with increasing film thickness. From the curvatures of the 
inverse susceptibilities it is thus possible to deduce the number of layers, which might be a way to extract information 
on the film thickness from experiment. For T < Tc, the inverse hard axis susceptibilites Xx~x are constants, but their 
value increases with increasing layer thickness, although not very strongly, which allows one to discriminate between 
films with differing numbers of layers. From Xxx(T — Tc), one can get in principle information about the exchange 
anisotropy strength, whereas the value of the Curie temperature depends on both the exchange interaction and the 
exchange anisotropy strengths. All the results discussed might be modified by a layer-dependence of the exchange 
interaction and exchange anisotropy or, when a different lattice type has to be considered. Other effects like domain 
wall motion or vortex excitations, which are not treated in the theory above, could also lead to modifications. 

The calculations here demonstrate that we are technically able to calculate the magnetic properties of in-plane 
anisotropic ferromagnetic multilayer films with S > 1/2. Hopefully, some of the predictions of the present calculations 
can be verified experimentally in the future, in particular with respect to the hard axis susceptibility. This might 
be possible if the experimental techniques discussed in RefA, where experimental results on a bilayer are reported, 
can be improved. High precision measurements of anisotropic susceptibilities of thin films, particularly above Tc, are 
called for, which is certainly a challenge for experimentalists. 

We are indebted to P.J. Jensen for useful discussions. 



* Also at Institut fur Theoretische Physik, Freie Universitat Berlin, Arnimallee 14, D-14195 Berlin, Germany.; Electronic 

address: froebrich@hmi.de 
' Electronic address: kuntz@hmi.de 

1 P.J. Jensen, S. Knappmann, W. Wulfhekel, H.P. Oepen, cond-mat/0303320 (2003), accepted for publication in Phys. Rev. B. 

2 P. Frobrich, P.J. Jensen, P.J. Kuntz, A. Ecker, Eur. Phys. J. B 18, 579 (2000). 

3 P. Frobrich, P.J. Kuntz, M. Saber, Ann. Phys. (Leipzig) 11, 387 (2002). 

4 P. Frobrich, P.J. Kuntz, Eur. Phys. J. B 32, 445 (2003). 

5 A. Ecker, P. Frobrich, P.J. Jensen, P.J. Kuntz, J. Phys.: Condens. Matter 11, 1557 (1999). 

6 P. Henelius, P. Frobrich, P.J. Kuntz, C. Timm, P.J. Jensen, Phys. Rev. B 66, 094407 (2002). 

7 W. Gasser, E. Heiner, and K. Elk, in 'Greensche Funktionen in der Festkorper- und Vielteilchenphysik', Wiley- VHC, Berlin, 
2001, Chapter 3.3. 



